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Abstract 

It is shown that an overcritical density plasma slab can be totally transparent to a p-polarized 
obliquely incident electromagnetic wave. High transparency is achieved due to the interference 
of the evanescent waves in the subcritical region. The transmission coefficient has the resonant 
character due to the excitation of a plasma surface mode (plasmon-polariton). It is crucial that two 
evanescent waves have certain phase shift to provide non-zero energy flux through non-transparent 
region. The required phase shift is obtained at the exact resonance and corresponds to the absolute 
transparency. In a warm plasma case, the excitation of the propagating longitudinal (electrostatic) 
modes becomes possible. We demonstrate that these longitudinal excitations facilitate the total 
transparency of an opaque plasma slab creating additional resonances in the transmission property 
of the system. 

PACS numbers: 52.27.-h, 52.25.0s, 52.38.Dx, 52.35.Lv, 52.38.-r 
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I. INTRODUCTION 



There exists a large number of theoretical and experimental studies dedicated to 
the problem of light interaction with dense media having negative dielectric permittiv- 
ity jll B B B ll^t • The increasing interest in the properties of such media has been driven by 
their potential applications in different branches of science and technology. Understanding 
of the fundamental processes that take place during the interaction of electromagnetic radi- 
ation with overcritical density plasmas (e = 1 — oj^/uj 2 < for uj < oj p ) is of great importance 
to a number of areas. In particular, one of the challenges of the inertial confinement fusion 
(ICF) experiments is the effective heating of the fuel pellet by the laser light (in the direct 
drive) or X-rays (in the indirect drive ignition experiments). The two well known absorp- 
tion mechanisms in plasma are of collisional P, Q| and resonant J, 0, ^| nature. In the 
collisional absorption process, the energy of light is converted into thermal electron motion 
through electron-ion collisions. In the resonant absorption mechanism, the linear conversion 
in nonuniform plasma leads to a resonant excitation of electrostatic plasma wave, which 
subsequently decays giving rise to energetic electrons. 

In this work we investigate another resonant mechanism that dramatically enhances the 
penetration of the electromagnetic field into overdense plasma. This mechanism involves 



surface plasma wave or plasmon-polariton 



' -I Q| which may be excited at the plasma- 



plasma interface. Unlike bulk plasma waves, which are electrostatic in nature, the surface 
waves are of mixed longitudinal-transverse type and can couple to radiating electromagnetic 
fields. As was shown in Ref . jl| , at certain conditions the energy of the incident electromag- 
netic radiation can be completely absorbed by the plasma through the resonant excitation of 
surface plasma waves which subsequently dissipate through standard irreversible processes 
like linear wave transformation or collisions. 

In this paper we show that the resonant excitation of surface modes may create the con- 
dition for absolute transparency of an overdense plasma slab. The transparency mechanism 
studied in our work is directly related to the problem of the anomalous transparency of 
light in materials with negative permittivities. Media with negative electric permittivity 
and magnetic permeability (^(J (the so-called left-handed materials or LHM) possess some 
unusual optical properties, in particular, the possibility of creating an "ideal lens" offering 
sub- wavelength resolution 0, Q, Q]- As was pointed out by Pendry jl7| |. the superlens 
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is essentially based on the possibility of evanescent waves amplification, facilitated by the 
excitation of surface plasmons. The amplification of the evanescent waves and their inter- 
ference are the basic mechanisms behind the anomalous transparency demonstrated in our 
work. This indicates the feasibility of superlensing in materials without any restrictions on 
the magne tic p ermeability fi, which has also been noted by Pendry. 



In Ref. 



2l| . the authors suggested that the incident electromagnetic wave may excite a 



pair of coupled surface plasma waves on both sides of the slab. The energy of the surface 
wave running on the opposite side of the slab can be dissipated by reemission of the elec- 
tromagnetic waves, which causes the layer to appear totally transparent. However a special 
symmetric structure consisting of media with different dielectric permittivities has to be 
formed in order for a pair of coupled surface waves to exist. 

In this paper we reexamine the conditions for total transparency of a layer of overcritical 
density plasma. Anomalously high transmission of light is explained through the interference 
of the evanescent fields formed by the incident /emitted electromagnetic radiation and the 
field of the surface plasma wave. We show that the total transparency of light can be 
achieved even when a single surface plasma mode is excited. Therefore, there is no need 
to form a symmetric structure in which an overcritical density plasma slab is sandwiched 
between two additional layers with specially chosen dielectric permittivities. In other words 
anomalous light transmission can be achieved by passing the light through a single plasma- 
plasma interface. In most laboratory experiments when high-power laser light interacts with 
a solid density target, a prepulse precedes the main laser pulse, creating the plasma density 
profiles discussed in this paper. 



II. SURFACE WAVE ASSISTED TOTAL LIGHT TRANSMISSION IN COLD 
PLASMAS 

As it was noted above, high transparency of dense plasma occurs as a result of the 
resonant excitation of a surface mode. Since the surface of the medium with e < behaves 
like surface wave resonator J2J, the problem of surface wave excitation is equivalent to that 
of resonator excitation. Under resonant conditions (external field frequency and its wave 
vector's tangential component have to match those determined from the dispersion relation 
for the plasma surface wave) the coupling between the surface eigenmode and the external 
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field becomes extremely efficient. The minimum in the reflection coefficient (the maximum 
in the transmission coefficient in the absence of dissipation) occurs for these resonant values 
of the external wave frequency and wave number. 

Here we review a dispersion relation for surface waves that exist on a cold plasma-plasma 
interface for arbitrary values of dielectric permittivities of the adjacent regions. It is found 
by solving the linearized set of cold fluid-Maxwell equations, 

V*E = -i§ ( la , 

f = " E ' < 1c > 
ot m 

where ions are assumed to be motionless. Letting all physical quantities vary as 
i&(x) exp (ik y y — iut), where x is the coordinate into the plasma and y runs along the 
interface with the field components of the form ~E=(E X , E y ,0), B=(0,0, B z ) one obtains, 

d / 1 dB z \ 2 ^ q 
dx \e dx J 

where a 2 = (k 2 — ^ej and e = 1 — uj 2 /uj 2 . Assuming that the interface is located at the 
plain x = 0, and for x < plasma has a constant density such that u p = u p i, and for x > 
a different constant density u p = uj p2) the solution to Eq. (J2J in both regions is, 

B z = Ai exp (aix), x < 
B z = A 2 exp (— a>2x) , x > 

It should be noted that here we only keep those solutions that decay away from the interface 
(a condition for existence of a well-defined surface mode). Using the boundary conditions 
at x = requiring the continuity of B z and dB z /dx/e, one readily obtains the dispersion 
relation for the surface waves on the cold plasma-plasma interface, 

aiAlW ^ ^ ^EE> ( 3) 

V 1 + A — 2ur 

where A = c^A/ciA, u is normalized to the plasma frequency in the region 2 where e < 
and k y is normalized to the classical skin depth 5 = c/uj P 2. As can be seen from Eq. (|5jL e\ 
and €2 must have opposite signs, so that one plasma region must have undercritical electron 
density while the other-overcritical. 



One important difference between the surface plasmons running on a plasma-vacuum in- 
terface (obtained from Q by setting A = 0)) and those on a plasma-plasma interface lies in 
a fact that the phase velocity in the latter case can exceed that of light, whereas the phase 
velocity of the former is always subluminal. Because of this fact, the surface plasmons on a 
plasma-plasma interface can couple to radiating electromagnetic fields, which are evanescent 
in plasma but may become propagating if a vacuum region is present beyond plasma lay- 
ers. The opposite is also true: if an external p-polarized electromagnetic pulse is obliquely 
incident on the structure comprising two plasma layers, it may under a resonant condition 
excite a surface plasma mode. On the contrary, a p-polarized electromagnetic pulse obliquely 
incident on a vacuum-plasma interface cannot excite a plasma surface wave since the phase 
velocity of this collective mode can never exceed the speed of light (the perturbation of the 
plasma surface by the external field propagates in this case with phase velocity greater than 
that of light) and the field of the external radiation is not evanescent (as we mentioned ear- 
lier, the fields of the surface wave are decaying away from the interface, or only evanescent 
external fields can excite the plasma surface mode). This means that if a layer of dielectric 
is added, so that a system comprising dielectric medium, vacuum and overcritical density 
plasma is formed in such a way that the external field in vacuum is evanescent (through 
total internal reflection of light in dielectric material) and the perturbation of the plasma 
surface by the evanescent field propagates with the phase velocity that is lower than the 
speed of light, then the external electromagnetic field may excite the surface plasma mode 
on a vacuum-plasma interface 0, |f| . 

Optical properties of a double-plasma system are determined by matching the fields of 
the incident /reflected electromagnetic waves with those of the surface wave. Let us consider 
a system which is characterized by spatial density (dielectric constant) distribution shown in 
Fig. ([Tjl. Ap -polarized plane electromagnetic wave is obliquely incident on this structure. In 
both plasma regions, the electromagnetic field B = (0, 0, B z ) is the solution to the following 
equation: 





with a general solution having the following from, 



B z = A ie ikx + A 2 e 



—ikx 



(5) 




The electromagnetic fields in the vacuum regions 
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{x < —d and x > a) have the following form, 

B z = E e ikocos8 x + R e - ikocos8 x , x<-d 



B z = Te ikocose x , x>a 



where T and R designate the field components for the transmitted and reflected waves. 
Matching solutions at different boundaries by requiring continuity of B z and dB z /dx/e 
across interfaces, one obtains a system of six equations (not shown here) for the unknown 
coefficients. The expression for the transmission coefficient has the following form, 

2 



T 



4e'(fcirf+fc2"-("+d)fcQ co5fl )A;ofcifc2eie2Cos0 

((l+e 2ifc l d )jfcoA:iteeie2 cos e+e 2ik 2 a (l+e 2ik i d )k k 1 k 2 e 1 e2 cos 6»-2ie i ( fe l d +* ; 2'»)x) 

T = k 2 €2 cos[/c 2 a] (kf + k^e{ cos 2 9) &m[k\d] + 



(6) 



{kiei cos[kxd\ (k\ + k%e\ cos 2 9) — ik [k\e\ + k\e\) cos 9 sin [kid]) sin[/c 2 a] 



where fci 2 = hoy/ (e^2 — sin 2 0). For resonant values of the external electromagnetic field 
frequency and incidence angle, one should expect anomalously high transmission of light. 

Figure (J2J) shows the transmission coefficient as a function of the incidence angle. As one 
can see, there is a sharp increase in the transmission properties of the system when the angle 
of incidence matches one of many resonant values of 9 = 0.671955 radians at which point 
the transmission coefficient reaches the value T = 0.99999. The resonant values exactly 
correspond to those given by an expression Q for the dispersion relation of the surface 
plasma wave. Thus, the anomalous transmission occurs even for the system consisting of an 
undercritical density plasma layer adjacent to an overcritical density plasma slab, so that 



there is no need to form a sandwich-like structure as argued in Refs. ( 3 5 S E|) i n which a 
second surface plasma wave has to be excited on the opposite side of the overcritical density 
layer to achieve the same effect. In other words, the anomalous light transmission can be 
achieved through the excitation of a single surface plasma wave. It is interesting to note 
that for a sandwich-like system, there are generally two peaks in the transmission coefficient 
at two different but close angles of incidence j^. When the second layer of undercritical 
plasma is removed, or the thickness of the overcritical density plasma slab becomes too large, 
the two solutions degenerate into one, which corresponds to the surface plasmon running on 
a single plasma-plasma interface. Analogous effect is seen in quantum mechanics where the 
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energy level of a particle in two-humped potential is split into two energy states that are 
close to each other Q]. As the potential is transformed into a single well, the two energy 
levels degenerate into a single energy state. 

As one can see from Figures (j2J) and (@J) the light transmission is a non-monotonous 
function of the widths of both plasma layers with a single maximum reached at certain 
correlated values. These values are determined by the interference condition between the 
two non-propagating (evanescent) fields. The evanescent fields are expressed through the 
hyperbolic functions (anomalous transmission occurs when both k\ and k 2 take on pure 
imaginary values). It can be easily seen that the total field in both plasma regions constitutes 
a sum of two exponential functions, one decaying with distance e~ x and the other growing 
e x (x points in the direction from left to right). It is a known fact that a single evanescent 
(decaying) wave carries no power. It can be easily shown that the x component of the 
Poynting vector for this wave is S x ~ Re[k] = (k is a purely imaginary spatial decay 
constant). However, when the total field is a superposition of decaying and growing modes, 
E,B ~ Ai exp[— ikx] + A 2 exp[ikx] {k is purely imaginary), the x component of the Poynting 
vector S x becomes, 



Since A1A2 — A 2 A* is purely imaginary (it can also be zero when the incident wave is 
evanescent and no dissipation is present), Re[k{A\A 2 — A 2 A\)] 7^ 0, signifying that the 
superposition of the decaying and growing modes carry a finite power. We shall call this 
superposition as the interference of the evanescent waves. The total light transmittance of 
an opaque plasma slab occurs when the Poynting vector inside the slab is equal to that 
of the incident radiation, so that Relk^AiA^ — A 2 Al)} = S x0 . Figures (jSJ) and © show 
the spatial distribution of the magnetic field modulus and a;-component of time-averaged 
Poynting vector S x = \Re[E y B* z ] for the case when surface wave is excited on a plasma- 
plasma interface. As one can see, the field amplitude grows dramatically in the first plasma 
layer (the magnitude of the magnetic field at the interface is more than fourteen times 
of that in the incident wave) to offset its exponential decay in the second to the level of 
the incident light. The x-component of the Poynting vector inside plasma regions is equal 
to that of the incident radiation-the energy of the field is perfectly transferred through the 
overcritical density plasma slab. For comparison, Fig. (J7J) shows the magnetic field amplitude 
for the case when there is no surface mode excited at the interface. Only exponentially 
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decaying solution to the field equations is observed. It should be noted here that light can 
be transmitted through an opaque high-density slab even without surface waves present at 
the boundary, if the slab is thin enough. The effect of the evanescent wave interference is 
at play in this case too. The difference lies in the actual values of the expansion coefficients 
A\ and A 2 in both cases. When there is no surface wave present, the coefficients are such 
that Re[k{A\A2 — A 2 Aiy\ < S x q, whereas if the surface wave is excited at the interface, 
the condition Re[k(A\A2 — A 2 A\)] = S x q can be fulfilled and the absolute transparency of 
optically opaque plasma layer is observed. 



III. SURFACE WAVES ON A WARM PLASMA-PLASMA INTERFACE 

The results obtained in the previous sections are applicable to the cold plasma case, 
limiting its applications to mostly solid-state plasmas. In plasma states created through the 
interaction of the high power electromagnetic radiation with solid structures, the electron 
temperature can reach values in excess of 1 keV. Therefore, it is important to understand how 
the finite temperature modifies the dispersion relation of the surface wave and subsequently 
the anomalous light transmission through the dense plasma. 

Kaw and McBridejl^] have used hydrodynamic description to investigate the finite tem- 
perature dispersion relation for electromagnetic surface waves propagating on a plasma- 
vacuum interface. We generalize their treatment for the case of surface waves on a warm 
plasma-plasma interface. The dispersion relation is found from solving the linearized set of 
warm fluid-Maxwell equations, 

VxE=— — 

c at 

47rn e„ 1 <9E 
V x B = v H — 

c c at 

V-B = 

V • E = -Airen 
dn „ 

^ + V.n v = 



dt m n 



where the subscript indicates equilibrium quantities; vth is electron thermal velocity; ions 
are assumed to be motionless. The above system of equations constitutes a closed set for 
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the unknown electron density, hydro dynamic velocity and electromagnetic fields. Just like 
in the cold case, it can be readily solved assuming that density and field perturbations are 
in the form ^f(x) exp[i(k y y — ut)}, where x is the coordinate into the plasma and y runs 
along the interface. This gives the following expressions for the electromagnetic fields and 
electron density perturbation in both plasma regions assuming that the interface is at x = 
and omitting the explicit dependence on u> and k y , 

„ _ r . AneAil - 0?) . , „ „ , . . . . 

E(i) = D(i) exp[aixj ^ ^ ' (*S e y + 7i e x) exp[7ixj, x < (7a) 

7i ~~ a i 

E( 2 ) = D (2 ) exp[-a 2 a;] - 47re ^( 1 — . (ik y e y - 7 2 e x ) exp[-7 2 x], x > (7b) 

n(i) = Aexp^x], x < (7c) 

n(2) — -B exp[— 7 2 x], x > (7d) 

where A, B, D are unknown coefficients, (3 = vth/c, e y , e x are the unit vectors in y and x 
directions respectively, 7 = (k y + -^p—\ , a = (k y + i ° p J° j ; subscripts 1, 2 designate 
both regions. As one can see the electromagnetic field in a bulk plasma is represented as a 
superposition of a transverse and longitudinal perturbations that are completely decoupled 
from each other. At this point we would like to bring reader's attention to the fact that 
just as in the cold plasma case, solution (JJJ) describes decaying branches for the transverse 
fields as well as for a longitudinal field in the overcritical plasma region (described by the 
decay factor 72 in our notations). However, the density perturbation as well as the asso- 
ciated longitudinal electric field in the undercritical plasma region are described not by an 
evanescent /decaying branch, but rather by a propagating oscillatory function (7! takes on 
purely imaginary value in this region). The presence of non-decaying oscillations in the 
longitudinal field signifies the existence of a bulk plasma wave in the undercritical plasma 
region. The existence of the bulk plasmon stems from the linear wave conversion process 
that occurs at the plasma boundary jl4j]. In the conversion process, the surface plasma wave 
gives rise to a bulk longitudinal plasma mode that propagates away from the interface into 
the undercritical density plasma. 

The dispersion relation for surface plasma waves on a warm plasma-plasma interface is 
found from the boundary conditions on the electromagnetic fields (tangential components 
of both fields e.q. E y and B z should be continuous across the interface). The two additional 
boundary conditions are obtained from the integration of the continuity equation and the 
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equation of motion over x coordinate in the limits x = ±5 and setting subsequently 5 = 0. 
Assuming that the density perturbation does not have any singularities at the interface (no 
surface charge present), the two additional boundary conditions have the following form, 

n n 
no 1 n 2 

n y|i = 7i yi2 (8) 

These are the continuity conditions for the electron density (normalized to its corresponding 
equilibrium value) and flux perturbations across the interface. In order to solve for all six 
unknowns in the system, the above four boundary conditions have to be also supplemented 
by the requirement that the divergence of the transverse part of the total electric field given 
by Eqs. (JTj) be equal to zero. Solving the system of six algebraic equations, one obtains 

[{U 2 - A) y/1 ~ U* + W y + {UJ 2 - 1) y/A - UJ 2 + k 2 ] 



(co 2 - A) + k 2 y + A {u 2 - 1) + k 2 \ = (A - l) 2 u 2 k 2 , (9) 

where A = u^/u^, lo is normalized to the plasma frequency in the region 2 where e < 
and k y is normalized to the classical skin depth 5 = c/u P 2. Eq. (jH)) is the general dispersion 
relation for surface plasma waves running on a warm homogeneous plasma-plasma interface. 
It is interesting to note that matching electromagnetic fields across the interface couples 
their longitudinal and transverse components, so that the resulting surface plasma mode 
is of mixed longitudinal-transverse type. Taking the limit A — > reduces this dispersion 
relation to that derived by Kaw [12] for the case of plasma-vacuum interface. The cold 
plasma dispersion relation is recovered from the limit /3 — > 0. Numerical solution of Eq. (jHJ) 
shows that for /3 1)2 < 0.02, which corresponds to electron temperatures 204 eV or lower, 
the cold plasma dispersion relation (jSJ) yields solutions that are close to those given by the 
general expression (JHJ), neglecting small imaginary parts in the eigenroots for a case of a 
warm plasma. The presence of imaginary values in the dispersion relation is directly linked 
to already mentioned conversion process. Because of the surface wave transformation into a 
bulk plasmon, its energy leaks out of the surface mode into a bulk plasma mode, which can 
be viewed as an effective dissipation of the surface plasmon, making the dispersion relation 
complex. It should be noted that the complex conjugate of the roots (with simultaneous 
reversal of the real parts of to — > — u and k y — > —k y ) are also solutions to the dispersion 
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relation (0), which is physically related to the opposite process in which the bulk plasmon 
converts into the surface wave. 

As plasma temperature gets higher, the eigenroots obtained from Eq. (JUJ) can significantly 
deviate from those given in the cold plasma limit as can be seen from Fig. (|SJ). This fact 
should be taken into consideration when investigating a possible surface wave induced plasma 
heating channel in ICF experiments where the initial temperatures (due to slab interaction 
with laser prepulse) can exceed 1 keV for ultra-intense laser pulses. 

IV. TEMPERATURE EFFECTS AND THE ANOMALOUS LIGHT TRANSMIS- 
SION 

Let us consider a structure similar to that of the cold plasma case and characterized by 
the spatial density (dielectric constant) distribution shown in Fig. (Q), but with non-zero 
electron temperature. A p -polarized plane electromagnetic wave is obliquely incident on this 
structure. In both plasma regions, the electromagnetic fields B = (0, 0, B z ), E = (E x , E y , 0) 
are the solution to the following system of equations: 



where k y = u/c sin 9 and (3 2 = 3v^ h /c 2 (the factor of 3 originates from the adiabatic 
equation of state) and e = 1 — uj 2 /uj 2 . The system of equations (fTUj) can be readily obtained 
from the initial set of Maxwell and warm fluid equations. It should be emphasized that 
Eqs. (Ill))) are only valid for not too high electron temperatures, so that the electron mean 
free path is smaller than the spatial inhomogeneity of the field inside the plasma. Otherwise, 
the hydrodynamic treatment is not adequate and one has to invoke kinetic description, as 
in a problem of anomalous skin effect P. \\\ . A general solution to the above system is 




(10b) 



(10a) 



(10c) 



B z = A ie ikx + A 2 e 



—ikx 



(11a) 



E x = —pt- (A ie ikx + A 2 e~ ikx ) 
k e 

E y = — (A ie ikx - A 2 e~ ikx ) 




(lib) 



(lie) 
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where k = u/c, k = k Q ^ (e — sin 2 9), k p = k y(e — (3 2 sin 2 9)/ (3. Coefficients A 1>2 repre- 
sent transverse part of the total field and Ci j2 represent its longitudinal/electrostatic part. 
The electromagnetic fields in the vacuum regions (x < —d and x > a) have the following 
form, 



B z = E e ikocos 


9 x _|_ ^j>g — ik() cos 9 


X 


x < —d 


E x = -E sm9e ikocosd x - 


Rsm9e- ikocosd 


X 


x < —d 


E y = E cos 9e ikocose x - 


Rcos9e- lkocose 


X 


x < —d 




-q Tg^ cos ^ 


X 


x > a 


E x = 


-Tsin9e lkocose 


X 


x > a 


Ey ~- 


= Tcos9e ik ° cosd 


X 


x > a 



where T and R designate the field components for the transmitted and reflected waves. By 
matching the solutions at all boundaries (requiring continuity of E y , B z , noV x and n/no 
across interfaces), one obtains a system of ten equations (not shown here) for the unknown 
coefficients. It should be noted that in Ref. [t| the authors used similar approach to inves- 
tigate the anomalous transmission of light through a structure consisting of an overdense 
warm plasma slab spaced between two underdense plasma layers. However, the choice of 
boundary conditions on the fields used by the authors, specifically the continuity of the spa- 



tial derivative of the x-component (z-component in their notations) of the electric field along 

2 

or reflection R 



E 



Ho 



x leads to non-physical solutions in which the transmission T 
coefficients can be greater than one. Boundary conditions employed in current work lead 
to solutions that always satisfy the identity |i?| 2 + |T| 2 = 1 (which is the conservation of 
energy for a plain- wave case for a non-dissipative system). The resulting expressions for the 
transmission/reflection coefficients are extremely bulky functions of electron temperature, 
electromagnetic wave incidence angle, its frequency, both plasma layers electron densities 
and thicknesses. 

As mentioned in the previous section, the dispersion relation for the surface plasma waves 
at non-zero temperatures can significantly deviate from that obtained in the cold plasma 
limit. Since the anomalous transmission of light is due to the excitation of the surface plas- 
mons, the resulting transmission coefficient exhibits thermal effects that are not seen in cold 
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plasmas. Figure Q shows the dependence of the transmission coefficient on the thermal 
electron velocity vth = \/T/m e . As one can see, it shows high frequency oscillations in the 
electron temperature space. The presence of these regular oscillations in the transmission 
of the system is related to the linear wave conversion processes that occur at the plasma 
boundaries. The surface wave excites a bulk longitudinal plasma mode propagating in the 
underdense plasma layer away from the plasma-plasma boundary. When it reaches the 
vacuum-plasma interface, the longitudinal plasma wave converts into a transverse propagat- 
ing electromagnetic wave, which interferes with the reflected light. Thus, in the general case 
the interference between the incident /reflected fields, surface plasma wave induced fields and 
the fields created through the linear conversion processes determine the transmission. 

The temperature oscillations can be understood using the following simplified model. The 
surface plasma wave excites a bulk plasmon propagating at angle 9' to the surface normal. 
This angle can be found from matching the y components of both modes' wave numbers to 
give: 

sin 0' = 4= sin (12) 
v e 

The plasma wave propagates the distance L = d/ cos9' with phase velocity v p h = c(3/\/e. 
Therefore, the phase difference A0 = uAt between the reflected light and that emitted 
through the conversion process is, 

A0= ! 1 =• (13) 
^/l-^sin 2 ^ 

When the phase difference is equal to a multiple of 2ir, both components interfere construc- 
tively giving rise to large reflection and thus low transmission. Therefore, one arrives at 
a following temperature condition at which minimal light transmission occurs through the 
system, 



P = \j ^ecsc* 9 - 6 CSC2 ^ Sln -, n=l,2,3,.... (14) 



With a simple estimate given by Eq. (fT4*j) . we are able to recover some of the minima in the 
transmission coefficient obtained from the complete solution of Eqs. 0. Other minima in 
the transmission coefficient, which are not described by Eq. (fT4*|) are due to other conversion 
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processes like the incident wave transformation into a longitudinal plasma mode at the 
vacuum-plasma interface that travels to the plasma-plasma interface, gets reflected (exciting 
a surface mode in a process) and travels back to the plasma-vacuum interface converting 
to the transverse electromagnetic wave. All these conversion processes contribute to the 
transmission properties of the double-plasma system, making it highly oscillatory function 
of not only the electron temperature but also the underdense plasma layer thickness d as 
can be seem from Fig. (|10|). The multiple peaks present in the transmission coefficient in 
this case can also be explained through the constructive interference between the reflected 
electromagnetic field and the field created as a result of the linear wave conversion process. 
Since the phase difference A0 also depends on the underdense layer thickness, one can arrive 
at a similar condition for the plasma thickness d at which minimal light transmission occurs 
through the system, 



A simple estimate given by Eq. (|15|1 can also recover the minima in the transmission coeffi- 
cient as a function of the undercritical density plasma layer thickness for a case of non-zero 
electron temperature. As the plasma temperature increases, the resonant condition for the 
excitation of the surface wave changes according to the dispersion relation given by Eq. Q, 
shifting the peak in the transmission coefficient to smaller angles. A distinct difference be- 
tween the zero temperature and a warm plasma case lies in a fact that there is only a single 
peak in the transmission coefficient (as a function of the undercritical density plasma layer 
thickness) in the cold plasma case, whereas the number of transmission peaks in a warm 
plasma is quite large. Unlike the cold plasma case where the energy of the external radiation 
is transferred through the constructive interference of evanescent fields, in a warm plasma 
case the superposition of the propagating longitudinal plasma wave fields (in a undercriti- 
cal density plasma slab) and the evanescent transverse electromagnetic fields determine the 
amount of the transmitted light as can be seen from figure where the x components of 
time-averaged total (which includes the contribution of plasma electrons) Poynting vector 
and the real part of the electric field are shown. As expected, the electromagnetic field has a 
form of the superposition of the evanescent and propagating/ oscillatory branches and under 
resonant condition the energy is totally transferred through the system. In addition, the 



d 




n = 1,2,3,.... 



(15) 
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results of our calculations show that in a warm underdense plasma layer the energy is mainly 
carried by the longitudinal electron plasma oscillations, facilitating the total transparency 
of the considered dual-layer system for quite large thicknesses of an undercritical density 
plasma slab. 



V. CONCLUSIONS 



It is shown here that overcritical density plasma can be made transparent to external 
radiation. The total transparency can be achieved even with a single surface mode excited 
on a plasma-plasma boundary. The anomalous transmission is explained through the in- 
terference effect between the evanescent fields of the incident electromagnetic radiation and 
those associated with the surface plasmons. Superposition of two (decaying and growing) 
evanescent modes can carry the power of the incident radiation over to the other side of a 
high-density plasma slab, making it totally transparent to light in a case when irreversible 
dissipation is absent. We have also shown that the excitation of the longitudinal plasma 
waves creates a multiplicity of additional resonances corresponding to the absolute trans- 
parency. We suggest that the mechanism investigated in this paper may be responsible for 
the anomalous transparency of metals observed recently in a number of experiments 8. Il7|. 

The role of dissipation effects is twofold, depending on possible applications of the surface 
wave phenomenon. In ICF experiments where high absorption of the incident radiation is 
sought, the presence of irreversible dissipation is necessary for effective heating of a fuel 
pellet through damping of the surface plasma waves. If dissipation is absent, the energy 
of the incident light will be totally re-emitted by these waves back into vacuum. If on the 
other hand a high transmission of light is needed, the dissipation plays a detrimental role, 
limiting the transmitting ability of the system. Therefore, a close control over a number of 
physical parameters (electron density, temperature, thickness) of the double-layer system is 
needed for successful technological implementation of the surface wave effect. 

The final point we would like to mention is related to the fact that since the Maxwell 
equations are invariant under reversal of time and space coordinates, the order at which 
the incident electromagnetic wave encounters both plasma layers (an undercritical density 
layer followed by that with overcritical density or vice versa) does not effect its transmission 
ability. 
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FIG. 1: Schematic diagram of the spatial density (dielectric constant) distribution. 
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FIG. 2: Transmission coefficient as a function of incidence angle. A double-layer plasma system 
becomes completely transparent at the incidence angle 9=0.671955. Undercritical and overcritical 
plasma slabs dielectric constants and thicknesses are e\ = 0.3428, d = 27 * c/oo P 2 and e2=-2.97, 
a = 3.12 * c/uj P 2 correspondingly. External radiation frequency oj/uj P 2 = 0.5019. 
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FIG. 3: Transmission coefficient as a function of a underdense plasma layer thickness. Overcritical 
plasma slab dielectric constant and thickness are £2=-2.97, a = 3.12*c/lo P 2 correspondingly. Under- 
dense plasma layer dielectric constant is e\ = 0.3428. External radiation frequency uj/uj P 2 = 0.5019. 
Incidence angle 6 = 0.671955. 
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FIG. 4: Transmission coefficient as a function of an overdense plasma layer thickness. Over- 
critical plasma slab dielectric constant is 62=-2.97. Underdense plasma layer dielectric constant 
and thickness are ei = 0.3428 and d = 50 * c/u P 2 correspondingly. External radiation frequency 
oj/uj p2 = 0.5019. Incidence angle 9 = 0.671955. 
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FIG. 5: Spatial distribution of an absolute value of the magnetic field associated with a single 
surface wave present on a plasma-plasma boundary. Undercritical and overcritical plasma slabs 
dielectric constants and thicknesses are t\ = 0.3428, d = 27 * c/lo P 2 and e2=-2.97, a = 3.12 * c/uj P 2 
correspondingly. External radiation frequency and angle of incidence are uj/lo P 2 = 0.5019, 8 = 
0.671955 correspondingly. Bq is the magnitude of the magnetic field in the incident wave. 
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FIG. 6: Spatial distribution of the x-component of the Poynting vector S x , normalized to the 
x-component of the Poynting vector of the incident wave S x q = 1/2-Eq cos^. 
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FIG. 7: Spatial distribution of an absolute value of the magnetic field when there is no surface 
wave present on a plasma-plasma boundary. Undercritical and overcritical plasma slabs dielectric 
constants and thicknesses are ei = 0.3428, d = 27 * c/lo P 2 and e2=-2.97, a = 3.12 * c/ui p 2 corre- 
spondingly. External radiation frequency and angle of incidence are u)/oj P 2 = 0.5019, 9 = 0.75 
correspondingly. Bq is the magnitude of the magnetic field in the incident wave. 
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FIG. 8: Dispersion relation for surface waves on plasma-plasma interface. The dashed line is the 
numerical solution of Eq. © for A = 0.3, f}\ = 02 = 0.1. The solid line corresponds to the cold 
plasma limit with A = 0.3. uj is normalized to the plasma frequency of an overcritical density 
region and k y is normalized to its classical skin depth c/iv P 2 
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FIG. 9: Transmission coefficient as a function of electron temperature (3 = ^j2>T/m e c 2 . Undercriti- 
cal and overcritical plasma slabs dielectric constants and thicknesses are ei = 0.3428, d = 27*c/uj P 2 
and €2=-2.97, a = 3.12*c/uj P 2 correspondingly. External radiation frequency and angle of incidence 
are oj/uj P 2 = 0.5019, 9 = 0.671955 correspondingly. 
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FIG. 10: Transmission coefficient as a function of a underdense plasma layer thickness d for plasma 
temperature parameter j3 = 0.055. Overcritical plasma slab dielectric constant and thickness 
are €2=-2.97, a = 3.12 * c/uj P 2 correspondingly. Underdense plasma layer dielectric constant is 
ei = 0.3428. External radiation frequency lo/u) P 2 = 0.5019. Incidence angle = 0.671955. 
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FIG. 11: Spatial distributions of the x-components of the Poynting vector S x (solid line) and 
the total electric field (dashed line) inside plasma layers for plasma temperature parameter (3 = 
0.15. Overcritical plasma slab dielectric constant and thickness are e2=-2.97, a = 3. * c/uj P 2 
correspondingly. Underdense plasma layer dielectric constant and thickness are e\ = 0.3428, d = 
37.5 * c/uj P 2 ■ External radiation frequency uj/uj P 2 = 0.5019. Incidence angle 6 = 0.64516. 
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